Noise Correlations in low-dimensional systems of ultracold atoms 
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We derive relations between standard order parameter correlations and the noise correlations in 
time of flight images, which are valid for systems with long range order as well as low dimensional 
systems with algebraic decay of correlations. Both Bosonic and Fermionic systems are considered. 
For one dimensional Fermi systems we show that the noise correlations are equally sensitive to spin, 
charge and pairing correlations and may be used to distinguish between fluctuations in the different 
channels. This is in contrast to linear response experiments, such as Bragg spectroscopy, which are 
only sensitive to fluctuations in the particle-hole channel (spin or charge). For Bosonic systems we 
find a sharp peak in the noise correlation at opposite momenta that signals pairing correlations in 
the depletion cloud. In a condensate with true long range order, this peak is a delta function and 
we can use Bogoliubov theory to study its temperature dependence. Interestingly we find that it 
is enhanced with temperature in the low temperature limit. In one dimensional condensates with 
only quasi-long range (i.e. power-law) order the peak in the noise correlations also broadens to a 
power-law singularity. 



I. INTRODUCTION 

The ability to trap ultracold atoms in tightly confined 
tubes formed by an optical lattice has opened the door 
for the controlled study of one dimensional physics. Such 
systems have been used to investigate ground state cor- 
relations [H, [il as well as dynamics [3j| and transport Q 
of strongly interacting bosons in one dimension. One 
dimensional traps of ultracold fcrmions have also been 
realized. The ability to control the interactions has 
been demonstrated using s-wave Feshbach resonances for 
fermions with spin and with p-wave resonances for 
spin-polarized fermions @ . 

From the theoretical viewpoint, one dimensional sys- 
tems provide good starting points to study strong correla- 
tion physics. Because of the enhanced quantum fluctua- 
tions, continuous symmetries cannot be broken in generic 
one dimensional systems, and mean field theories fail. 
Nevertheless the Luttinger liquid framework provides a 
well developed formalism to treat these systems theoret- 
ically (see for example Q)- In place of a mean field order 
parameter the tendencies to ordering manifest themselves 
by slowly decaying algebraic correlations and correspond- 
ingly divergent susceptibilities. A diverging susceptibility 
in a particular channel implies that coupling an array of 
tubes in the transverse direction would lead to true order 
in that channel. Thus weakly coupled one dimensional 
systems provide a possible theoretical route to investigate 
open questions regarding competing orders in higher di- 
mensions. 

But the absence of long range order, which makes 
one dimensional systems interesting also complicates the 
ways by which these systems can be probed. Various 
methods have been proposed and used to probe specific 
correlations. For example time of flight [l| as well as in- 



terference experiments [8| can probe single particle cor- 
relations. Bragg d, and lattice modulation Q spec- 
troscopies can measure dynamic density correlations . In 
[lll | we proposed that noise correlations can be used as a 
highly flexible probe of ID Fermi systems. In particular 
this method is sensitive to a wide range of correlations 
including spin, charge and pairing, and it treats the vari- 
ous channels on an equal footing. These results, obtained 
using the effective Luttinger liquid theory were confirmed 
by Luscher et. al. using a numerical simulation of a mi- 
croscopic model. 

In this paper we provide the detailed theory of quan- 
tum noise interferometry in one dimensional Fermi sys- 
tems and extend it to interacting Bose systems. 

Atomic shot noise in time of flight imaging was pro- 
posed in Ref. [13] as a probe of many body correlations 
in systems of ultra cold atoms. Specifically what is mea- 
sured using this approach is the momentum space corre- 
lation function of the atoms in the trap 

Gaa'{k,k') = (riakna'k') - ("-afe) ("a'fc' ) 7 (1) 

where a, a' are spin indices, fc, k' are momenta, and Ua^k 
and Ua'^k' are the occupation operators of the corre- 
sponding momentum and internal state. This relies on 
the assumption that the atoms are approximately non- 
interacting from the time they are released from the trap 
(even if they were rather strongly interacting in the trap). 
Here we envision a system of one dimensional tubes, that 
allows expansion only in the axial direction. To ensure 
weak interactions during time of flight, the radial con- 
finement in the tubes can be reduced simultaneously with 
release of the atoms from the global trap. 

Recent experiments demonstrated the ability to detect 
many-body correlations by analysis of the noise correla- 
tions ll]). For example long ranged density- wave correla- 
tions (induced by external lattice) were seen in a system 
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of ultracold bosons [l^flU and fermions[l^ in deep opti- 
cal lattices. The sign of the correlation, peak or dip, de- 
pended on the statistics as expected, demonstrating that 
the effect is essentially a generalization of the Hanbury- 
Brown Twiss effect. Pairing correlations in fermions orig- 
inating from a dissociated molecular condensate were also 
observed in experiment [l7|. All of these demonstrations 
involved states with true long range order that are easily 
related to simple two particle effects such as Hanbury- 
Brown Twiss or bound state formation. However, one of 
the main points of Ref. [3] was that the noise corre- 
lations can be used to measure more general correlation 
functions in many-body systems, even in the absence of 
true long range order (LRO). Here we demonstrate this 
idea by developing a detailed theory for the noise correla- 
tions in one dimensional systems of bosons and fermions. 

The analysis is laid out in the following order. In sec- 
tion HJ we consider a model system of spinless fermions 
in one dimension. To elucidate the connection between 
noise correlation and ordering tendencies we first assume 
in section lll Al the presence of true long range order. This 
allows to provide a clear physical picture for the way or- 
der in spin, charge (particle-hole) or pairing (particle- 
particle) channels translate into sharp features (in this 
case delta function peaks) in the noise correlations. In 
section III Bl we shall consider the actual system of in- 
teracting fermions in one dimension where only power 
law order parameter correlations exist (quasi-long range 
order). This system is characterized by a tendency to 
charge density wave (CDW) ordering for repulsive inter- 
actions, and to superconducting order (SC) for attrac- 
tive interactions. The naive expectation is that the delta 
function peaks will be replaced by power-law singularities 
with a power that reflect the algebraic order parameter 
correlations. Interestingly we find that such a simple 
relation exists only if there are order parameter correla- 
tions that decay with a sufficiently slow power. This is 
the case beyond a critical interaction strength (repulsive 
or attractive). At weak interactions, on the other hand, 
the situation is more subtle and singular signatures of 
both the dominant and sub-dominant ordering tenden- 
cies are seen in the noise correlations (see Fig. [2l3,c). We 
compare these results with the information that can be 
extracted from measuring the static structure factors in 
the spin and charge channels. 

In section Hm we move on to analyze the noise correla- 
tions in one dimensional systems of interacting spin-1/2 
fermions. In the long wavelength effective field theory the 
interaction is parameterized by two Luttinger parameters 
Kp and K„ corresponding to the spin and charge sectors, 
as well as a backscattering parameter gi_L. The ordering 
tendencies of this system are summarized in the phase- 
diagram in Fig. [3] (See [11]). We find that the noise 
correlations provide direct information on the real space 
order parameter correlations only for negative backscat- 
tering. In this case a spin-gap opens and the system is 
characterized by competing CDW and SSC correlations, 
which are revealed by the noise correlations. 



In section IIVI we turn to interacting Bose systems. 
We first study the noise correlation within Bogoliubov 
theory for a condensate with true long range order (sec- 
tion |TVX|. The interesting feature her is a peak in the 
noise correlations ^ at k -\- k' = Q. This can be simply 
understood as the signature of pairing correlations in the 
Bogoliubov wave- function, describing quantum depletion 
of the condensate. We compute the evolution of this sig- 
nature with increasing temperature. In addition we point 
out the existance of sharp dips in the noise correlations 
along the lines fc' = and fc = 0, these reflect correla- 
tions between the condensate and the quantum depletion 
cloud. In section HV Bl we move on to discuss interacting 
one dimensional Bose liquids, where the single particle 
density matrix decays as a power-law with distance. We 
show how the sharp features in the noise correlation func- 
tion on the lines fc + fc' = 0, fc — 0, and fc' = broaden 
to power-law singularities. 



II. SPINLESS FERMIONS 

A system of interacting spinless fermions is perhaps 
the simplest Fermi system with a non trivial competi- 
tion of ordering tendencies. Such systems can be im- 
plemented in experiments using fully polarized fermionic 
atoms that interact via an odd angular momentum chan- 
nel. The interaction strength can be tuned by using a 
p-wave Feshbach resonance @. Another possibility is 
to use a Bose-Fermi mixture, where the phonons of the 
bosonic superfluid mediate effective interactions between 
the fermions [l^ . 

In the restricted geometry of ID the Fermi surface con- 
sists only of two points, the left (L) Fermi point at —kf 
and the right (R) Fermi point at -t-fc/. We introduce the 
left- and right-moving fields ^fj^ and ip^i^ through: 

^(x) = e-'^^^^VL + e^'^^V-i? (2) 

i)L and il}R are slowly varying fields, because the rapidly 
oscillating phase factors e^'*^^^ have been separated out. 
The natural order parameters that characterize this sys- 
tem are the charge density and superconducting (pairing) 
operators: 

OcDW = i^R^^L = P2kf (3) 

Osc = i>Ri^L (4) 

The central problem we address in this section concerns 
the connection between the correlations of these physi- 
cal order parameters and the observable noise correlation 
signal. 

A. Long range order 

The relation between order parameter correlations and 
sharp features in the noise correlation function is most 



3 



apparent when the system supports true long range or- 
der. Of course one dimensional Fermi systems generically 
do not display spontaneous long range order in either the 
CDW or SC channels. However long range order can be 
induced by an external potential (e.g. a superlattice po- 
tential) or it can form spontaneously in a weakly coupled 
array of one dimensional systems. Moreover the intuition 
gained from the exercise will be useful in approaching the 
more interesting case of power-law order parameter corre- 



lations (quasi long range order), which will be considered 
in the following sections. 

Consider first the case of long range pairing correla- 
tions (i.e. SC order). Take k (k') to be near the right 
(left) Fermi points, and q = k ~ kp , q' = k' + kp devia- 
tions from those Fermi points. It is now useful to write 
the noise correlations explicitly in terms of the pairing 
correlations 



{rikrik') ~ J2 I ci^ti^'drdr'e'(«+«')(^-^')e'(«~«')(''-'^''/2(v4(X + r/2)V'l(X-r/2)^i(X'-r72)V^fl(X' + r72)) 



-^2 j dXdX'drdr'e*(''+«')(^"^')e'(«-'')('^-''')/2^Bt(X)B^-(X')) 



(5) 



r 



Here the operator Br{X) creates a fermion pair whose 
constituents are separated by r and their center of mass 
coordinate is X. Note that we have dropped terms that 
undergo rapid 2k p spatial oscillations and therefore van- 
ish under integration. The existence of long range pairing 
correlations implies 



(6) 



at long distances (\X — X'\ oo). Here $(r) — 
(■!/'■'■ (X r/2)-0t(X - r/2)) is the translationally invari- 
ant pairing wave function. The long range saturation of 
the correlation function leads to a singular contribution 
to (O at g = -q': 



{nku-k') ^ S{q + q )\^q)\^' 



(7) 



We thus conclude that the noise correlations in this case 
are directly related to the long distance limit of the pair- 
ing correlations. This is one way to formally justify the 
mean-field decoupling of the noise correlation function 
carried out in Ref. which gives: 



(riknk') - (nkUk') = |(^ki^_k) |'<5(k + k') (8) 



A very similar analysis follows for the case of long range 
order in the particle-hole cannel. Take for example CDW 
order at the wave- vector 2kp. In this case we should write 
the noise correlation function in terms of the density cor- 
relations to expose their singular contribution: 



{ukUk') = --^ f dXdX'drdr'e*(«-«')(^-^')e*(«+«')('-^''')/2(Vj]j(X-|-r/2)^L(X-r/2)i/it(^/„^//2)^^(X' + 



= --!/' dXdX'drdr'e*(''+'')(^-^')e'(''-^')('--^')/2(p2;..,,.(X)p2fe.,.'(X')) 



(9) 



Note the minus sign in front of the integral, which re- 
sulted from commuting the Fermi operators to obtain an 
expression written in terms of density wave correlations. 
Long range CDW correlations imply 

{p2k,AX)p2k,,r'{X')) ^ S(r)S(r'). (10) 

which leads to a singular contribution to ^ for q ^ q' 

{nkUk') -^-S{q-q')\E{q)f (11) 

This justifies the corresponding mean field decoupling of 



the noise correlation function: 



inkHk' 



lit 



t 

k+2kp 



tk)\'S{k-k') 



(12) 



Again, as in the SC case, we see that the noise correla- 
tions are directly related to the long distance limit of the 
non decaying order parameter correlation. 

We note two obvious distinctions between the cases of 
order in the particle-particle (pairing) and particle-hole 
channels. In the former the singular correlations are be- 
tween particles with opposite momenta and are positive 
correlations. In the latter case, by contrast, the singular 
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FIG. 1: Schematic representation of several types of fluctua- 
tions. In these diagrams, k and k' represent the momenta of 
the atoms and holes, and q and q' the momenta relative to 
the Fermi points, (a) A pairing fluctuation (or Cooper pair), 
which is the dominant fluctuation in the SC phase. These fluc- 
tuations result in positive noise correlations along q — ~q' . 
(b) A particle-hole (p-h) fluctuation associated with a CDW 
state. This fluctuation results in negative correlations along 
q — q' . (c) Two p-h pairs. This fluctuation results in both 
positive and negative correlations. Positive correlations for 
q = ~q' and negative for q = q' . 



correlations are between particles whose momenta dif- 
fer by 2kF and are negative, that is anti-correlations. To 
better understand the origin of this effect it is worthwhile 
to inspect the mean field wave functions, which sustain 
the respective broken symmetries. For the superconduct- 
ing order the mean field wave-function is the BCS state 
representing a pair condensate: 



■^sc) ^ll[uk + vkipiijlk\ |0) 



(13) 



Contrary to a filled Fermi sea the particle number Uk 
in a specific k point is not definite in this wave function. 
However, if a particle is found at fc, there is with certainty 
another one at —k. This implies positive correlation be- 
tween Tifc and as visualized in Fig. [T] (a). 

The CDW state on the other hand may be viewed as 
a condensate of particle-hole pairs on top of the filled 
Fermi sea. It is then written in a way that exposes the 
similarity to the BCS state: 



■2k, 



\FS] 



(14) 



Here I FS ) denotes a full fcrmi sea wavefunction. As 
in the BCS state, the particle number Uk at a specific k 
point is not definite. Now however, if one finds a particle 
at k > kp there will be with certainty a missing particle 
(hole) in the Fermi sea at fc — 2fci?. This implies anti- 
correlation between Uk and nk-2kF as visualized in Fig. 
[11(b). 



B. Quasi-order 

Because of strong quantum fluctuations, an actual one 
dimensional Fermi system cannot sustain true long range 
order such as the SC and CDW orders discussed above. 
Instead, at zero temperature a critical phase with power 
law correlations, or quasi long range order, is established. 
If the power-law decay is sufficiently slow then it is rea- 
sonable to expect that it would still make a singular con- 
tribution to the integrals in ([5]) or ([9]). To calculate the 
resulting singularity and its dependence on the system 
parameters we use the low energy Luttinger liquid the- 
ory. The basic idea is that the asymptotic low energy 
and long wave-length properties of the interacting one 
dimensional Fermi system are captured by a universal 
harmonic theory: 



S = 



— [ dxdT\c(d^ef +c-^(dre?] 
tK J ^ 

J dxdr [c(a^$)2 + c-\dr^)^] (15) 



2ttK 
K 
2^ 



Here K is the Luttinger parameter, which determines 
the power-law decay of long range correlations, c is a 
sound velocity, which we will henceforth set to 1. The 
bosonization identity, which relates the bosonic fields to 
the fermion operators is: 



1 



\/27ra 

with the commutation relation: 



[eix),m] = Jiog^. 

2 a — IX 



(16) 



(17) 



a is an artificial cutoff of the bosonization procedure 
which must be sent to zero at the end of the calculation. 
Implicit in the action ([15]) is also a physical short distance 
cutoff xq, taken in most cases to be of order l/kp- Here 
$ is the phase field of the Fermi operator. H = is 
the smooth (fc « 0) component of the density fluctuation, 
and by (fTT]) . is also canonically conjugate to $. 

The effective action (fTS]) is a free theory in terms of the 
bosonic fields. It therefore allows to calculate the needed 
correlation functions using simple gaussian quadrature. 
The fermion interactions affect the correlation functions 
only through the Luttinger parameter K. For non inter- 
acting fermions K — I, K < 1 for repulsive interactions 
and K > 1 for attractive interactions. In general K de- 
viates more from K = 1 the stronger the interactions. 
However it is in general not possible to make a precise 
connection between the microscopic parameters and the 
Luttinger parameter. Among other things, our analysis 
points to a way of extracting this parameter from exper- 
iments. 

For any value of K the system shows either CDW 
or SC quasi-long range order, which means that the 
correlation functions decay slow enough to give a di- 
vergent susceptibility. The calculation of these cor- 
relation functions has been given in numerous places 
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(see e.g. 0, HH). For the long distance behav- 
ior of OcDW we have {OcDwix,T)OcDwiO,0)) ~ 
cos(2fc/)(a;2 + c'^r^y{^-acnw)/2 ^ ^^^^ fp^. q^^ ^^^^ 

{O sc [x, t)Osc (0,0)) - (a;2 + c2r2)-(2-"sc)/2 ^^cre the 
scahng exponents for CDW and SC are 



OiCDW — 2 

asc = 2 



(18) 
(19) 



The susceptibilities correspond to the spatial and tem- 
poral Fourier transform of the correlation functions, 
X{k,uj) ~ / dxdTe'''^+"^^ {O{x,t)O{0,0)). These will be 
divergent at large distances exactly if the scaling expo- 
nent of the operator 0{x, r) is positive. As we can see 
from (fT8|) and (fT9|) . xcdw diverges at /c = 2fc/ for if < 1, 
and xsc diverges at fc = for if > 1. In this sense of 
QLRO we say that the system is in the CDW reg ime for 
if < 1, and in the pairing regime for if > 1 [21[, as de- 
picted in Fig. [2]e). We note that from Eq. [T5l and [TBI 
one can read off the duality mapping: 0<-><I>, if<-^if~^, 
which leaves the action invariant, and maps the CDW 
regime < A' < 1 onto the SC regime, 1 < if < oo. 

With this formalism we derive the noise correlations 
G{q, q'), which we show in detail in App. |Al We find 



Giq,q') 
with 

and 



piqxi2+iq' X34 

{27t)^L •^(^i2)-^*(^34)(.4-l)(20) 



a — IX 



A - 



' {Xq + x14){Xq + X23) 

- (a^o + ^13) (2^0 + 2^24) 



(21) 



(22) 



The integration in Eq. [^D] is over the three spatial vari- 
ables a;i2, X23 and 2:34. The exponents g and h are given 
by g = (if + if-i - 2)/4 and h={K- if"^)/4. 

Following the discussion of the MFA, we note that this 
integral 'contains' the equal-time correlations of Ocdw 
and Osc- This can be seen by setting a;i4 = and 
2^23 = 0, which gives the Fourier transform in q — q' of 
-{OcDw{^)OcDw{'^))- If we set a;i3 = and a;24 = 0, 
we obtain the Fourier transform of {O sc{^)0 sci'^)) in 
q + q' ■ To discuss this further we introduce the following 
variables: z — {xi2 — a::34)/2, — (xn + X23)/'2, and 
/i_ — {xii — X23)/2. With this, [uqUqi), is of the form 



pi{q-q' )z+i{q+q' 

T{z + h^)T*{z - h^) 

A^ + (/^+ + /^_)^A2 + (/^+-fe_)^ '» 
A^ + {z + h+y K^ + iz-h+Y ) ' 



In the limit if <C 1, we have /i ^ 0. This enforces the 
integrand to be negligible except in regions with z, /i_ w 
0. With this, the expression approximately evaluates to 
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FIG. 2: a) - d) Noise correlations {nqUgi) — {nq){nqi) of a ID 
Fermi system for different values of the Luttinger parameter 
K. For a) - d): if = 0.4, 0.8, 1.25 and 2.5. For if < 
1/2, we find a negative algebraic divergence along q — q' , 
reflecting a quasi-condensate of p-h pairs. For if > 2, we 
find a positive algebraic divergence along q — —q', indicating 
a quasi-condensate of p-p pairs. In the intermediate regime 
we see signatures of both p-h and p-p correlations, which are 
singular at the origin, e) 'Phase diagram' of a Luttinger liquid 
of spinless fermions as defined by diverging susceptibilities 
and signature of the order parameter correlations in the noise. 
For strong interaction, if < 1/2 and if > 2, we find diverging 
noise correlation along a line. For if < 1/2 it is a negative 
peak on the line q — q' signaling a particle-hole correlation. 
For if > 2 it is a positive peak at g = —q' signaling a particle- 
particle correlation. In the intermediate regime (1/2 < if < 
2) there are signatures of both the leading and sub-leading 
fluctuations in the noise correlations which diverge at g = 
q' = 0. 



In the dual limit, if ^ 1, we have h ^ 0, which en- 
forces z, « 0. In this limit the integral approximately 
evaluates to 



sgn(2if'^ 



l)\q + q- 



i\2K- 



(25) 



{uqUq,) - sgn(2if-l)|g-g'| 



/|2i<'-l 



(24) 



which can also be inferred from duality. These contribu- 
tions are divergent for if < 1/2 and if > 2, and turn out 
to be the dominant contribution in theses regimes. 

To confirm this expectation, and in order to study the 
regime 1/2 < if < 2, we evaluate this integral numeri- 
cally for different values of if. In Fig. [5] we show Q{q, q') 
for K = 0.4, 0.8, 1.25 and 2.5. The K = 0.4 example 
shows indeed a power-law divergence of the particle-hole 
type, which we find througout the if < 1/2 regime. For 
if > 2 we find a result similar to the if = 2.5 exam- 
ple, an algebraic divergence of the particle-particle type. 
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These two regimes are indicated in Fig. [2]e). 

In between these two regimes, for 1/2 < K < 2, we 
find a regime in which both p-h as well as p-p correlations 
exist, i.e. we find precursors of the near-by competing or- 
der. A simple argument for the qualitative shape of the 
noise correlations function in this regime is the follow- 
ing: If we consider an interacting ID Fermi gas and treat 
the interaction perturbatively, the lowest order contribu- 
tion would consist of states that contain two p-h pairs, 
as two fermions have been taken from the Fermi sea and 
put into the unoccupied states above the Fermi sea (see 
Fig. [l]c)). Such a state exhibits qualitatively the noise 
correlations that are observed in Fig. [2]b) and c): posi- 
tive correlations for q < 0,q' > and q > 0,q' < 0, and 
negative correlations for g > 0, g' > and q < 0,q' < 0. 
To quantify how this affects the line-shape of the noise 
correlations we expand G{q, q') to second order in h (see 
App. |B]for details). The result is 

g{q,q') ~ -hhgn{q)sgn{q')mm(^,j^). (26) 

This expression shows a divergence at g — q' ~ 0, 
and both particle-particle and particle-hole fluctuations. 
Higher order terms either enhance particle-particle (for 
K > 1) 01 particle- hole {K < 1) fluctuations. 

Before concluding the section, we compare the static 
and dynamic structure factor with the noise correlations 
discussed in this section. The dynamic structure factor 
can be measured by a stimulated two-photon process, as 
has been demonstrated in [9]. It is defined as 

S{k,u;) ^ {p{k,u:)p{^k,-u)). (27) 

As discussed in Ref. [Tj] , it is given by 

S{k,uj) - K\k\5{Lu - c\k\) (28) 

for small fc « 0. For k k, 2nkf, it behaves as 

S{2nkf + q,uj) - [u^ -c^q^y'^-^ (29) 

where n ^ 0. So at ii" = 1, coming from larger values, 
two additional peaks appear at ±2A:/, because the ex- 
ponent K - 1 in Eq. ([^ switches sign. At K = 1/4 
two additional peaks appear at ±4fc/, and so on. So 
the measurement of the dynamic structure factor, being 
the Fourier transform of the density-density correlation 
function, certainly allows insight into the CDW phase, 
but contains no information about the SC phase. In 
contrast, as demonstrated in Fig. [2l the measurement 
of noise correlations proposed in this paper allows the 
identification of pairing and CDW ordering in a single 
approach and on equal footing, because both the CDW 
and the SC correlations are contained in the expression 
for the noise correlations. 

The static structure factor of a BEC has been mea- 
sured in [To] . It corresponds to the instantaneous density- 
density correlations: 

S{k) ^ (p^kPk) (30) 



For fc « 0, we have 

S'(fc«0) - K\k\ (31) 

and for k « 2nkf: 

S{2nkf + q) ^ |g|2"'^-i. (32) 

Here the first set of peaks appears at if = 1/2, the next 
set of peaks at K = 1/8 and so on. This is the power- 
law divergence that dominates the noise correlations for 
K < 1/2, as discussed in the previous section. Again 
there is no signature of SC. 

In summary we have derived the noise correlations of a 
spinless fermionic LL in this section, and compared it to 
a MFA result. We found different subregimes with qual- 
itatively different behavior in each of the quasi-phases, 
summarized in Fig. [2l and discussed how phenomena 
such as QLRO and competing phases are reflected in the 
noise correlations. 



III. SPIN-1/2 FERMIONS 

In this section we discuss the noise correlations of an 
SU(2)-symmetric Fermi mixture (20| . 

Our analysis applies to systems of the form 

H = -tY, ^l,>j.- + C^E"^.T"U (33) 

{ij).a i,(J 

i.e. the ID Hubbard model, or a mixture in a ID contin- 
uum with a contact interaction: 

^ = E/ + ^ / "T(^)^i(^) (34) 

which can both be realized in experiment. 



A. Long range order 

Before we turn to the LL picture, we introduce the 
types of order that occur in this system, and discuss what 
kind of signature can be expected if they develop long- 
range order, with similar arguments as in Sect. Ill Al 

As for the spinless fermions we split the spinful 
fermionic operators into left- and right-moving fields: 

V'T/i(x) = e-^^/^VL,T/i+e^'^"V'i?.T/i- (35) 
We introduce these order parameters 0, : 



Oss = ^RA'^L.x (36) 

Ots = ■ipR.^'^LA (37) 

OsDW = V'ljjV'La (38) 

OcDW = V^lijV'L,!- (39) 
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Oss describes singlet pairing, and therefore contains t 
and I operators, whereas Ots describes the x and y com- 
ponent of triplet pairing, and therefore describes pairing 
between equal spin states. Osdw is the order parameter 
of the X and y component of the spin density wave order. 

Because these order parameters contain both R and 
L, as well as "f and | operators, we expect correlations 
between and ?iL,f, as well as nn^i and n^,!, so we 
consider two types of correlation functions; correlations 
between atoms in the same spin state 

Gr^{q,q') = ("t,9"t,?'> - K,9>K,9'> (40) 

and correlations between opposite spins 

Qui<l,Q') = {nLqn^,q') - {ni,q){n^^q,). (41) 

We again use the same convention that q is located near 
the right Fermi point and q' located near the left Fermi 
point. By considering the dominant contraction of these 
correlation functions, analogous to Eq. ([8]) and Eq. p^ . 
we can expect the following signatures: For CDW order 
we expect {q, q') ^ Sq^qi , for SDW Qi^ {q, q') ^ -5q,q' 
whereas for triplet pairing we expect G^'[{,q,q') ~ ^q,-q'^ 
and for singlet pairing Giiiqjq') ~ ^q,-q'- Each of these 
statements can be confirmed and quantified with a mean- 
field calculation. 



B. Quasi-order 

We bosonize these fermionic fields according to: 



1 



= ±i©T/-l(2:)+i*T/i(2;) 



(42) 



with the same definitions of Oyi and that we used 
for spinless fermions. We introduce spin and charge fields 
according to: 



:(eT±ex) 



(43) 
(44) 



Written in terms of and dp{x) — Qp{x) — kfX the 
action of any SU(2)-symmetric system separates into a 
charge and a spin sector: 



S — Sp + 5*0- 



with: 



and: 



S n 



2ttK, 



Vp{dc 



py^x'Jpj 



1 j ]_{^drQaY +V„{d,Q,f 



cos(V8^e<,) 



(45) 
(46) 

(47) 
(48) 



1, 



Singularity at origin 



p-h/p-p 


p-p/p-h 


(SDW/CDW) 


(TS/SS) 


1/2 


1 2 


CDW 


CDW/SS-cusp 


SS/CDW-cusp 


] 

ss 


(CDW) 


(CDW/SS) 


(SS/CDW) 


(SS) 



FIG. 3: "Phase Diagram" summarizing the singular signa- 
tures in the noise of the various order parameter correlations 
in a one dimensional spinfuU Fermi system. gi± is a backscat- 
tering parameter and Kp the luttinger parameter in the charge 
sector. The divergent susceptibilities in the different regimes 
are in brackets. For positive backscattering only weak sig- 
natures of order parameter correlations are seen through a 
weak singularity at the origin. On the other hand in the spin 
gapped phase at negative backscattering the noise correla- 
tions give detailed information on both the CDW and SSC 
correlations. 



Each of these sectors is characterized by a velocity Vp/^- 
and a Luttinger parameter Kp/^j. In addition to the 
quadratic terms in the action we find a non-linear term, 
describing backscattering processes in the spin sector, 
with the pref actor gi^±. If this action is derived from a 
system with a short-ranged interaction such as p3p and 
P4p . these parameters have the following properties: For 
repulsive interaction, one finds Kp < 1, K„ > 1 and 
positive backscattering 51, j_, for attractive interaction 
Kp > 1, Ka < 1, and negative backscattering gi^± < 0. 
As discussed in [2l|, j^^, this sine-Gordon model, can be 
treated with an RG calculation to identify two limiting 
cases: The case in which the backscattering term is irrel- 
evant and the system flows towards the non-interacting 
fixed point (K^- ^1), which happens for repulsive inter- 
action, and the case in which the backscattering term is 
relevant and a spin gap appears (K^ — *■ 0), which hap- 
pens for attractive interaction. In the evaluation of the 
noise correlation functions we will use these limiting val- 
ues, Kcr — 1 for the gaplcss phase, K^ — for the spin- 
gapped regime. One can find the phase diagram of this 
system in exact analogy to the spinless case, by studying 
the correlation functions of the order parameters. The 
scaling exponents of these operators are given by: 



ass 
ars 

O^SDW 
OLCDW 



K„ 



2 - K^' - K, 
2-Kp- K-' 
2-Kp-K„. 



(49) 
(50) 
(51) 
(52) 



From these expressions one can read off the structure of 
the phase diagram. In the gapless phase we have K^ = 1 , 
therefore singlet and triplet pairing, as well as SDW and 
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FIG. 4: Noise correlations {n-^^qTi^ qi) — {ni^q){n-j ,ji) (a - 

c), W,9"-T,9'> - {"-i,9>("'T,9') (d - f). and {ntot,qntot,q') - 
{ntot,q){ntot,q') is - i) of a spin-1/2 Fermi system in ID in 
the spin-gapped phase for different values of the Luttinger 
parameter Kp. Kp — 0.8 for (a, d, g), Kp = 1 for (b, e, h), 
and Kp = 1.25 for (c, f, i). In (a-c) we can see CDW order- 
ing in the t|-channel, in (d-f) we see singlet pairing in the 
It-channel. The noise correlations in the total density clearly 
show the coexistence of orders. For Kp < 1 CDW is domi- 
nant, and singlet pairing is subdominant, for Kp > 1 it is the 
other way around. 

CDW are algebraically degenerate. For Kp > 1 we find 
a TS/SS phase, for Kp < 1 we obtain a SDW/CDW 
phase. For Kg- 0, both asDW and ars are sent to 
— oo, whereas ass and acDW are now given by ass — 
2 — K~^ and acDW = 2 — Kp. Hence, we can distinguish 
four regimes: For Kp > 2 we have singlet pairing, for 
Kp < 1/2 we get CDW ordering. In between these two 
values of Kp the system shows coexisting orders, that 
is, both the singlet pairing susceptibility and the CDW 
susceptibility are divergent. For 1/2 < Kp < 1 CDW is 
dominant and SS is subdominant, for 1 < Kp < 2 it is 
the other way around. 

The noise correlations can be calculated in the same 
way as described for the spinless case. We obtain anal- 
ogous expressions to Eq. (j20p . in which the exponents g 
and h are replaced by: 

5TT/iT = (Kp + K;'+K, + K;' -4)/8 (53) 

and: 

h^^ = {Kp-K;')/8+{K,-K;')/8 (54) 

/iXT - {Kp-K;')/8-iK,-K-')/8 (55) 

In order to understand in what regimes of the phase di- 
agram we should expect algebraic divergencies, we again 
consider the equal-time correlation functions of the op- 
erators ([55)1 - ([55)1 . In momentum space, these correlation 
functions scale as where a is the corresponding 

scaling exponent, given in ([49|) - ((52l) . If the system is 
in the gapless phase (i.e. K^^ — 1), these correlation 
functions never exhibit an algebraic divergence, and we 
should expect to find coexisting fluctuations throughout 



the phase diagram for gi^± > 0, similar to the regime 
1/2 < K < 2 for spinless fermions. If the system is in 
the spin-gapped phase, we find the following behavior: 
Both TS and SDW fluctuations are frozen out, i.e. only 
short-ranged, whereas SS and CDW are increased by 1, 
compared to the gapless phase, because K^ — > 0. We 
therefore expect algebraic divergencies for Kp < 1 in the 
It channel, and for Kp > 1 in the ti channel. 

A numerical study of the noise correlations confirms 
these expectations: We indeed find coexisting fluctua- 
tions in the gapless phase for any value of Kp. Further- 
more, since the expressions (j54p and (j55|) become identi- 
cal for Ka- = 1, we find that G^iiq,q') = Giiiq,q') in this 
regime. This is a manifestation of the degeneracy (at the 
algebraic level) of triplet and singlet pairing for Kp > 1, 
and of spin density and charge density wave ordering for 
Kp < 1, as discussed in ^21]. For the spin-gapped phase 
(K^ — > 0) this symmetry is broken and and be- 
have qualitatively different. We find the following be- 
havior: G^iq,q') shows an algebraic divergence of the 
p-h type for Kp < 1, as expected from the equal-time 
correlation function of the CDW order parameter, an al- 
gebraic cusp for 1 < Kp < 2, and no ordering for Kp > 2. 
Ql-\{q,q') behaves in a complementary way: an algebraic 
divergence of the p-p type is found for Kp > 1, a,n alge- 
braic cusp for 1/2 < Kp < 1 and no ordering below that. 
In particular, for Kp in the vicinity of 1, we find coexist- 
ing orders, as we demonstrate in Fig. [3) This is particu- 
larly clear if we consider the noise correlations of the total 
density ntot.q = n^,q + ni^q, for which the noise correla- 
tions are given by Gtotiq,q') = 2Q^^{q,q') + 2Qi^{q,q'). 
In Fig. dig) - i) we clearly see the coexistence of pairing 
and CDW ordering. 

To understand this behavior further in this limit, we 
use the same argument as for the K < 1/2 and K > 2 
regimes for spinless fermions. We re- write {n-\ ,qn-\ ,qi) and 
{ni^qTii^q') in the same way as (|23p . where g and h need 
to be replaced by giy-^i and h-^y^^, respectively. In the 
limit Ka- the arguments that lead to the scaling 
behavior ([M]) and (pS)) become exact, and we obtain 

{n^,qn^.,q') ~ sgn(ifp-l)|g-g'|^'--i (56) 
(n^^qn^.q,) ^ Sgn(if;l - l)\q + q'\''^"-\ (57) 

These expressions show exactly the structure that was 
found numerically: algebraic divergencies for Kp < 1 
{Kp > 1) in the ff (t|) channel, and an algebraic cusp 
foil<Kp<2 (1/2 <Kp< 1). 

IV. BOSONS 

We turn to address the noise correlations in bosonic 
systems with either long range or quasi long-range or- 
der in the off diagonal density matrix {b^{x)b{0)). As 
in the fermion case we shall start with the case of true 
long range order, relevant to three dimensional systems 
at temperature T < Tc. This analysis is also relevant for 
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lower dimensional systems if they are sufficiently weakly 
interacting. Then the fact that the off diagonal density 
matrix decays as a power law is unnoticeable in a con- 
densate of realistic size. We shall discuss in some de- 
tail the possibility of observing the pairing correlations 
associated with quantum depletion and show how such 
measurements would depend on the temperature. Then 
we move on to derive the noise correlations in a one di- 
mensional Bose system at zero temperature, taking into 
account the power-law behavior of the correlations. As in 
the case of Fermions we use the effective Luttinger liquid 
theory, which correctly accounts for the singular contri- 
butions to the noise correlations due to the long distance 
power-law behavior of the off diagonal density matrix. 




0.01 0.1 1 10 



FIG. 5: Noise correlation of pairs G{k, — fc) of a BEG, as a 
function of k^, for different ratios of At/C- 



A. Bose-Einstein condensate with true ODLRO 

Our starting point for analysis of the noise correlations 
is the Hamiltonian of a weakly interacting Bose gas with 
contact interactions 

H = 51 ^^ka^ak + ^ X! "k+q«P-q"p«k- (58) 

k 

Here ek = k'^ /2ra is the free particle dispersion with m 
and U — Anas/m with the s-wave scattering length. 
To compute the correlations in the condensed phase we 
apply the standard Bogoliubov theory (see e.g. [1^). As 
usual, the operators aj and oq are replaced by a number 
representing the condensate amplitude ^/Nq, while the 
other modes are treated as fluctuations and expanded to 
quadratic order. The effective Bogoliubov Hamiltonian 
is then given by 

Hb = J2 + Upo) a[ak + ^ E «-k«k' (59) 

k k 

where po = Nq/V is the condensate density. This Hamil- 
tonian is diagonalized by the Bogoliubov transformation 
flk = UkCkk + fka^k with = (cjk + Sk + ?7po)/2c^k, 
vl = (-Wk + ek + U po)/2wk, and Wk = V<^k(ek + 2Upo). 

The structure of the ground state wave function in the 
Bogoliubov approximation is given by 

I *B ) - exp y^al + ^(wk/uk)aLk4 I ) • 

(60) 

Like the BCS state ([T3l) . the Bogoliubov wave- function 
describes perfectly correlated pairs of particles at mo- 
menta k and — k, which suggests the appearance of pair- 
ing correlations in the noise. 

It is straight forward to compute the noise correlations 
G(k, k') for k, k' ^ 0. Because the Bogoliubov Hamil- 
tonian ([55)1 is quadratic we can use Wick's theorem to 
decouple the four point function 

^(k,k') = (nk)(l -I- (nk))^kk' + |(aLka|.>l'^k,-k' (61) 



where the expectation values correspond to thermal aver- 
ages, and Tik — ajj.ak. A bit more care is needed if either 
k = or k' = because the quadratic hamiltonian de- 
scribes only the fluctuations in the depletion cloud, not 
in the condensate number. To obtain the fluctuations 
in the condensate within Bogoliubov theory, we use the 
conservation of total particle number, which implies that 
fluctuations in the condensate number are exactly to mi- 
nus those of the depletion cloud. In other words we sub- 
stitute no = N — X^k^i^o ""k f*^^ the condensate particle 
number operator. Then we may use (|59|) to compute the 
noise correlation between points fc = and k' . Putting 
it all together we get the general expression for the noise 
correlations: 

Q{k, k') = gk<5k,-k' + fkSkk' - hi^iSko + ^k'o) 

+'5k0^k'0 

q 

where 

5k = ulvl (1 + 2{na,k)f 

/k = .9k + ('^a.k) (1 + (?^a,k)) 

hk = 25k + ("-Q,k) (1 + ("-a,k)) (63) 

Here (riQ,k) — (o^kak) = [exp(wk/r) — 1]^"'^ is the quasi- 
particle number distribution. 

At zero temperature each term in (|62p has a simple 
physical interpretation. We already noted that the first 
term manifests the pairing correlations present in the 
quantum depletion described by the Bogoliubov wave 
function (j60p . The second term is a positive correlation 
due to boson bunching at a point in k-space. The dips 
at 5(k, 0) reflect the fact that an extra atom found at 
k in the quantum depletion cloud corresponds to a pair 
of atoms, now missing from the condensate. Finally the 
positive peak at C/(0, 0) appears because extra atoms in 
the condensate must always come in pairs, annihilated 
from the depletion cloud. That is, if we find an extra 
atom in the condensate, then we are sure to find another 
extra atom in it. 
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The evolution of the peaks with temperature and their 
momentum dependencies are controlled by the ratio of 
two natural length scales of the problem. One is the 
healing length of the condensate which is determined by 
interactions S, = l/y/2mUpo. The other is the thermal 
wavelength Xt — 1/ \/2mT. The combination u^v^^ is of 
course temperature independent and may be expressed 
using the healing length alone as [8{^k)'^ + 4:{^k)'^]~^. 
The temperature dependence arises from the ratio Wk/T, 
which appears in the distribution function, and may be 
expressed as {Xt / 0^ (k^) \/ 1 + (fc^)^. Note that the di- 
mensionless ratio is the same as the ratio n/T. 

An interesting correlation to observe is the pairing cor- 
relation at k' = — k, which can be written explicitly using 
(PI and dMl) as 



5(k,-k) = 



coth^ 



wr)(A;e)yr+W? 



sm^+ii^kr 



(64) 



This function is plotted on a log-log scale in Fig. [S] 
The behavior of the peak weight at small relative mo- 
mentum fc^ << 1 has a very simple form. First, at 
very low temperatures, such that {fj,/T)k^ » 1 we 
have tj(k, — k) w l/[8(fc^)^]. This is the power law seen 
for T = in Fig. [5l At higher temperature or suffi- 
ciently small momentum such that {fi/T)k^ << 1 we 
have ^(k,-k) w (T/^)V[2(fc^)'*], which is seen in the 
other curves in the same figure. 

It is interesting to note that the pairing correlations 
are substantially enhanced with temperature. This seems 
surprising given that the origin of the pairing is the quan- 
tum depletion in the ground state ([SO)) . The effect may 
be interpreted as Bose enhancement of paired thermal 
fluctuations. However we should also note that the over- 
all noise level is also growing with temperature, that is 
the local (in k) particle number fluctuation G'(k, k) is in- 
creasing even more steeply with T. In Fig. [S]we plot the 
pair correlation normalized by the local number fluctua- 
tion 



P(k) 



g(k,-k) 



(65) 



Clearly, for T^O (i.e. Xt/^ oo) we have P{k) = 1 Be- 
cause in the ground state ([SO)) the number fluctuations 
always come in opposite momentum pairs the particle 
number at k and —k must be identical, and hence also 
the correlations G{k, —k) and G{k, k) are equal. At finite 
temperature the pairing correlations are suppressed com- 
pared to the local (in k) number fluctuation. The ratio 
remains 1 at very high momentum because at Wk » T 
thermal occupation (which is exponentially suppressed is 
negligible compared to the quantum depletion. More in- 
teresting is the fact that P(k) approaches 1 also in the 
limit of small relative momenta, which may again be a 
signature of Bose enhancement of the pairing fluctua- 
tions. 



P(k)^ 
0.75 

0.5 
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FIG. 6; Pair correlation function P{k) = Q{k, —k)/Q{k, k), as 
function of k^, for different ratios of At/C- 



B. Quasi-condensate 

In this section we discuss the noise correlations for a 
LL of bosons. We use Haldane's representation [2^ of a 
bosonic operator, defined as: 



h{x) = [p( 



'0 



(66) 



The fields Q{x), ^{x) and n(a;) are defined in the same 
way as for fermionic LLs. po is the average density. Note 
that now the sum is over the even harmonics, 2m, and 
not the odd ones, 2m -I- 1, that are used to represent 
a fermionic operator. The action of the system can be 
written as: 

S 5— / (fxdjd^'e - — / d?xda^d''^ (67) 

2ttK J 27r 7 

As for fermionic systems, these representations are just 
quadratic in the fields, therefore all correlation functions 
reduce to Gaussian integrals. 

For small momenta k and k' the noise correlations are 



G{k,k') 



pI 



i'''^''+'''^''T{xi2)T{x^^){Ah - im) 



which we derive in detail in App. |C| T and A are of sim- 
ilar form as before, with g = 1/{AK) and h = — \/{AK). 

This expression can be numerically evaluated for a fi- 
nite system, with the replacement x L sin(27rx/L)/27r. 
In Fig. [7] we show G{k, k'), plotted for different values of 
K. 

For large values of K, the LL results resemble quali- 
tatively the BEC result: The noise correlation function 
shows both a sharp pairing and bunching contribution 
which are equal in magnitude, as well as sharp quasi- 
condensate contributions along (fc,0) and {0,k'). In ad- 
dition there is a large peak at (fc = 0,fc' — 0), (invisible 
in the representation in Fig. [7]). As we reduce K to 
smaller values (corresponding to larger repulsive interac- 
tions between the bosons), the quasi-condensate gets vis- 
ibly broadened. This is due to the fact that there is no 
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Ok Ok 

FIG. 7: Noise correlations Q{k,k') of a ID Bose system for 
K = 1.05, 2.5, 8 and 20 (a-d). The grey scale is linear, with 
light colors corresponding to positive values, and dark colors 
indicating negative values. The center of the plot corresponds 
to (0, 0). For large values of K, the properties of a ID bosonic 
superfluid approach the ones of a BEG with various sharply 
peaked structures. For smaller values of K, the pairing peak 
and the quasi-condensate dips broaden and become weaker. 



Finally, we discuss the transition of a ID bosonic su- 
perfluid to a Mott insulator, which can occur if there is 
a lattice potential present that is commensurate to the 
density of the superfluid. For fc. A;' w and ior K ^ 0, A 
approaches 1, as can be seen from (1691) . In this limiting 
case, the integral of x becomes S^k — S^^k'- The remain- 
der of the integral can be evaluated to be (nfc)^ for k = k', 
so the entire singular contribution is given by {nk)'^Sk,k' ■ 
Therefore, for iiT — > 0, which describes the Mott insula- 
tor transition, we find Q{k,k') — > (nk)'^dkk' — {nk){n'k'), 
as for the higher dimensional case [Tj]. By using the 
higher modes of Haldane's representation, we can deter- 
mine the behavior of Q(k,k') for k « 2nkb and k' « 
2mki,. We use the expression b ^ y^e^™®e**, and 
we get: {uk) = Po J dx i2e'^'' J^n{x 12) , with g given by 
g = 1/AK + v?K. For (nfen^/) we obtain {ukrik') — 
plj e'''''^^+''''''^^Tn{xi2)T^{x3,i)A. A is defined as be- 
fore, with h now given hy h = —nmK — l/AK. We can 
now go through the same steps that were used to identify 
the (5-function for fc, k' « 0. We find that for any n and 
m, the noise correlation function approaches the peaked 
shape that was found for fc. A:' w in the limit K ^ 0. 
This coincides with the result found in [l^, in which an 
ansatz of the form \MI) = Y[i^l\0) was used to derive 
this result. 



V. CONCLUSIONS 



true condensate in ID: The occupation in fc-space is not 
peaked as a (5-function, but only an algebraic divergence, 
which gets broader for smaller values oi K. As a con- 
sequence, also the pairing peak gets broadened, because 
the distribution of the total momentum of the pairs that 
are created from the condensate gets broadened. Fur- 
thermore, the overlap with the broad quasi-condensate 
dips diminishes the magnitude of the pairing peaks for 
smaller values of K, and eventually, for K — 1, which is 
the Tonks- Girardeau limit, the pairing peak is entirely 
surpressed. 

In contrast to the broadening of the pairing peak and 
the quasi-condensate dips, the bunching peak is a S- 
function for all values of K. This arises because we 
have \g\ = \h\, and therefore the integrand in Eq. (jC5|) 
does not fall off in one direction. This can be under- 
stood from the integral expression (|C5|) . We introduce 
the variables x = {xi2 — X34)/2, z — (x\2 + 2;34)/2 and 
y = a^23 + -Z, and rearrange the integral (|C4[) : {ukUk') = 
/ Q^^k^+^k^ j:{^z + y)T{z — y)A. Here we introduced the 
definitions Afc ^ k — k' , k = k + k' , and 

X ^ / {xl + {x + z)^)ixl + {x ^ z)^) \h 

\{xl + {x + yr){xl + {x-y)^)J ' ^ > 

From this, it is clear that for 2: — s- 00 the integrand ap- 
proaches 1, and does not fall off to zero, giving rise to a 
(5-function. 



We have investigated the nature of noise correlations in 
systems of ultra cold atoms, that do not necessarily sup- 
port long range order in any order parameter. The essen- 
tial difference between noise correlations fl]) , which may 
be measured in time of flight experiments [14 [Tsl. Kdi. IVq. 
and standard order parameter correlations, is that the 
noise correlations are non local in real space. In other 
words, they are not a fourier transform of a two point 
correlation function in real space. In this paper we inves- 
tigated, the general connections, that nevertheless exist 
between these two types of correlations in many-body 
systems. We then derived the singular contributions to 
the noise correlation function for several ultra-cold bo- 
son and fermion systems of interest. Our focus was on 
low dimensional systems, which display power law order- 
parameter correlations. 

For Fermi systems, we showed in sec lira , that 
true long range order in a particular order parameter 
(spin, density or pairing) leads to a delta-function contri- 
bution to the noise correlations. It is tempting to assume 
that if there are power-law decaying correlations in the 
same order parameter, then they would contribute cor- 
responding algebraic singularities to the noise. However, 
we showed that this is only the case provided the alge- 
braic decay is sufficiently slow. For example in a system 
of one dimensional spinless Fermions, the noise correla- 
tions appear similar to the mean field result, albeit with 
power-law peaks, only for K < 1 (strong repulsive in- 
teraction) or if > 2 (strong attractive interaction). In 
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the intermediate regime 2 > K > 1/2, they display a 
unique structure that reveals both the dominant and sub- 
dominant order parameter correlations. 

In a system of spinfuU fermions, on the other hand, 
the noise correlations display significant singular contri- 
butions only in the spin-gapped phases (i.e. for negative 
backscattering) . In this regime they provide information 
on both the singlet-pairing and density wave correlations. 

In all cases we showed, that the noise correlations treat 
the particle-particle (pairing) and particle-hole (spin or 
density) channels on the same footing. Noise correlations 
are thus sensitive, and can easily distinguish between the 
dual order parameters. In general, pairing correlations 
manifest as positive noise correlations (i.e. peaks), while 
spin and charge correlations appear as anti-correlations 
in the noise (dips). This is in marked contrast to external 
probes, such as Bragg scattering, which couple only to 
spin or particle densities (particle-hole channel). 

When considering Bose systems in section HVl we first 
treated a condensate with true off diagonal long range 
order using Bogoliubov theory. The most interesting fea- 
ture in this case is a correlation peak ~ (5(k-|-k'), which is 
a direct manifestation of the pairing correlations in the 
quantum depletion cloud. We find, somewhat counter 
intuitively, that these correlations are enhanced with in- 
creasing temperature at low temperatures. It should be 
noted however that the " normalized" pairing correlation, 
that is, relative to the number fluctuation on the fc-point, 
is indeed suppressed with temperature. Additional peaks 
appear due to correlations between the condensate num- 
ber and the number of particles in the depletion cloud. 
These appear on the lines fc = and fc' = 0. 

Finally we used the quantum hydrodynamic (Luttinger 
liquid) description [2^ of the Bose liquid, to address one 
dimensional Bose systems at T = 0. For this system, 
characterized by algebraic decay of off diagonal order, 
the pairing correlation in the noise also converts to a 
power-law singularity. 
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APPENDIX A 

In this appendix we derive Eq. I20i using bosonization. 
We write {uqUqi) and (n,) as 

•' i=i 

(V'lj(l)^i^(2)V'i(3)^L(4)) (Al) 



(n,) = dld2e'''--(V'Jj(l)^H(2)). (A2) 

q is the momentum relative to the right Fermi point, 
q — k — kp, and q' the one relative to the left Fermi 
point, q' = k' + hp, as before. We use the notation T' 
for 'xi\ etc., and Xi2 for xi — X2, etc. An analogous 
expression holds for (nqi), with 'R' replaced by 'L'. 

Using the bosonization expression in Eq. (jl6p for (rig) 
we obtain 



rfld2 
27raZ' 



Jqxi2 /p-ieil)-i'S>{l) ie{2)+i-S>{2) 



{A3) 



and a similar expression for n^/ with 8(2) — > —0(2) and 
-9(1) ^ 6(1). By using e^e^ = gA+Sg[A,B]/2 gg^. 



dld2 
2'KaL' 



g»g2;i2/gi(e(2)-e(l))+i(*(2)-'I>(l))\ 



,[e(i),*(2)]/2+[*(i),e(2)]/2 



(A4) 



The commutator between the fields Q{x) and (f>(x) is in 
Eq. [171 Next we use (e'^) = e^^ Here it is necessary 
to impose a short distance cut-off xq on the interaction, 
that is, K needs to depend on the momentum and has 
to fall off exponentially to 1 for momenta of the order of 
I/xq. The expression that we use is 



((e(x)-e(o))^ 



1 



+ x^ 



K-l^ xl + x"^ . 



An analogous expression holds for $(a:) with K replaced 
hy K^"^ . The cut-off xo on the interactions will stay finite 
throughout the calculation and can be interpreted as an 
effective bandwidth. If we apply this to 'T-g(g'), we obtain 



{nq(q-) 

where we introduced 



= ^ yrfa;i2e*«^-.F(*)(xi2) (A5) 



with g given hy g ^ {K + K ^ - 2)/4. 



(A6) 



We can evaluate 



along the same lines. We 



again use Eq. 1161 rearrange the exponents in the same 
way as we did for Uq, while keeping track of the non- 
vanishing commutators between them, and take the ex- 
pectation value, to obtain 

„iqxi2+iq' xa 

^ -((e(i)-e(2)-e(3)+e(4))^)/2 
(27ra)2L2 

-((*(l)-*(2)+*(3)-*(4))^)/2 



g[e(i),$(2)]/2-H[$(i),e(2)]/2 

g-[e(3),$(4)]/2-[<[.(3),e(4)]/2 

This can be evaluated to 

r,iqxi2+iq' X3i 



(A7) 



{nqUq') 



(27r)2L 



-T{xi2)T*{x:ii)A (A8) 



13 



The integration in this expression is over the three spatial 
variables X12, X23 and X34. A is defined as 



A 



i4 + 4i)i4 + 43) \'' 

{Xq + Xi^){Xq + X24) 



(A9) 



The exponent h is given by h = {K~K ^)/4. Combining 
the expressions that we derived for (uq) and {UqUqi) we 
get for g{q,q'): 



giq,q') = 
which is Eq. \M 



(27r)2L 



-^(a:i2)^*(a:34)(^-l) 



APPENDIX B 

In this appendix we expand Q{q,q') to second order 
in the exponent h. We show that the first order term 
vanishes, and the second order term gives Eq. The 
first order term is given by 



h j dx2z log(^ 



[xl + S{xiif){xl + S{x23f) 



s{x 



13 



s{x2iy 



(Bl) 



with s{x) defined as s{x) = L/27rsin(27rx/L). Here we 
explicitly kept the expression for a finite-size system. 
Since the integral separates into the sum of four terms 
of the form ± J'log(x§ + s(a;)2), and each of these four 
terms integrates to the same value (as can be seen by 
shifting the integration variable, and using the period- 
icity of sin(27ra:/L)), the entire expression integrates to 
zero. The second order term is given by 



2 



j dx23 log ( 



' {Xq + xl^){xl + X23) 

• i^o ~^ ^13) (^0 ~^ ^24) 



(B2) 



where we left out the cut-off and the finite-size represen- 
tation for notational convenience. It can be checked that 
these will ensure the following expressions and manipu- 
lations to be well-defined. 

Expression (jB2[) corresponds to a sum of integrals of 
two types: / dx\og^{x^) and /(r) = J dx log x^ \og{x + 
r)2, as can be seen by expanding the square in (|B2|) . The 
first integral merely provides terms that cancel divergent 
terms of the second integral type in the limit L ^ 00. 
The second term is given by /(r) = Ci + C2\r\, with 
Ci and C2 some constants. This can easily be shown by 
taking derivatives, and by observing that f{r) is even. 

With this expression for the various integrals of the 
type f{r), that we get from expanding the square in (|B2p . 
we get for g{q, q'): 



S{q,q') - 



aiqxi2 piq X34 



X12 X34 

-2|a;i2|-2|a;34|) 



-i\xi2 + X34I -I- 1X12 - 2^34! 

(B3) 



To evaluate this integral we divide the integration range 
into eight sectors. As an example we treat the case given 
by: X12 > 0, 2:34 > 0, X12 — X34 > 0. (The other sectors 
are characterized by similar sets of inequalities.) For this 
part of the integration range we get an expression of the 
form: 



piqxi2 1 

dxi2 — (e''"i^-l) (B4) 

X12 iq' 



The other sectors of the integration range give similar 
expressions. These can be grouped into a sum of inte- 
grals that contain integrals of the form J dx/x exp(ipa;) = 
«7rsgn(p), with different combinations of q and q' in the 
exponent. With that we get for G{q,q') an expression 
which can be written as 



G{q,q') 



which is Eq. 



-h'^sgn{q)sgn{q') min (^■^^, , (B5) 



APPENDIX C 



We calculate (rik) for fc « mode, for which the Bose 
operator is given by 6 ~ ^/poe**. For (rife) we find: 



{nk) ~ Po 



da;i2e^'=^-e-^((*(2)-*(i))^> (ci) 



For the evaluation of the expectation value ((<I'(2) — 
<I>(1))2) we use a slightly different cut-off procedure than 
for the fermions, in particular: 

(($(2)-$(i))2) = _Liog^i±^. (C2) 

With that we find (uk) ^ Po J da;i2e*'^^" J^(xi2). where 
we defined: 



(C3) 



The exponent g is given hy g — l/AK. Next we evaluate 
the expectation value {ukUk') along the same lines. We 
obtain: 



{ukUk') ^ Pq 



2 / ^ikxi2-\-ikx34 



H^12)H^34)A (C4) 



A is defined in the same as in Eq. (|A9[) . but h is now 
given hy h — —1/AK. We combine these expressions to 
get the correlation function G{k, k'): 



g{k,k') ^ plj e 



2 / ikxi2+ikx34 



T{X12)J'{X34){A-1()C5) 
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